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Abstract 

Commonly used techniques to study non-perturbative aspects of the strong interactions have a deep connection with 
rational approximants, and in particular with Pade approximants to meromorphic functions. However, only recently 
this connection has been acknowledged and efforts at fully exploiting it are only starting. In this article I will briefly 
review the most prominent techniques used in non-perturbative strong interactions with special emphasis on its relation 
with Pade approximants. I will then concentrate on a set of open problems outside the scope of these conventional 
techniques where Pade approximants might be extremely useful. 
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1. Introduction 



Particle physics is nowadays a mature discipline inside theoretical physics, where all observed phenomena among 
its elementary constituents (quarks and leptons) can be described in terms of four fundamental interactions: gravi- 
tational, strong, weak and electromagnetic. Unless one goes to extremely high energies, where gravitational effects 
are no longer negligible, strong, weak and electromagnetic interactions suffice to give a satisfactory description of the 
dynamics of elementary particles. Even though the nature of those three interactions is apparently very different, there 
is a common guiding principle, the gauge principle, that asserts that particle dynamics are determined by symmetries. 

In this article I will concentrate on the strong interactions, which describe the interactions of quarks mediated by 
gluons. Free quarks are described by the Dirac Lagrangian: 

£ = ifr[iy^-m]i/ f , (1) 

while interactions are implemented by requiring invariance of the previous Lagrangian with respect to local S U(3) 
transformations. This prompts the appearance of the gluon field and the final Lagrangian is 

1 6 

-Cqcd = - -Tr CTG^ + V $ 9 [iy.lT - m q ]if, q , (2) 

9=1 

where ZX 1 = & 1 — igsG 11 and [D M ,D V ] = -igsG^y, g s being the strong coupling constant. The kinetic term is added to 
make the gluon a dynamical field. Thus symmetry fully determines the structure of the strong interactions, with the 
6 quark species (flavours) belonging to the fundamental representation of S U(3) and the gauge potential to the 
adjoint representation. 
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The simplicity of Eq. (f2| does however little justice to the complexities it hides. First of all, quarks and gluons have 
never been observed. Instead, experimentalists only detect hadrons, and a large number of them. Therefore, Eq. (f2]i 
was initially received with scepticism and only taken as a model where hadrons were composite objects and quarks 
and gluons their building blocks. While this complied very well with some static properties of hadrons, there was no 
dynamical explanation of this compositeness or confinement: the equations of motion cannot be solved analytically, 
and perturbative expansions in the strong coupling g s seemed sheer nonsense. 

A major breakthrough came with the realization that at extremely high energies fi the binding of quarks and gluons 
becomes loose, i.e., 

lim g s Qx) = . (3) 

This asymptotic freedom I22I I41I1 of the strong interactions means that QCD at high energies behaves like a non- 
Abelian version of quantum electrodynamics (QED), and in particular makes perturbative expansions in the strong 
coupling (pQCD) meaningful. 

Parallel to the development of pQCD, techniques were developed to grasp some of the non-perturbative effects of 
QCD. The operator product expansion (OPE) Il50i 14711 is a well-defined procedure to incorporate non-perturbative 
effects as inverse powers of momenta q 2 . In practice, given a quantity Tl(q 2 , g s ), at large values of q 2 one has a double 
expansion of the form 

m\ 8s) = J] g'l log"' q 2 L + 2 2g] , (4) 

n.m V jfc " ) 

where each of the series is believed to be asymptotic. 

At very low energies, there is also (non-perturbative) information that can be extracted. The so-called chiral sym- 
metry is broken and Goldstone theorem requires the presence of massless particles 

EL 

whose interactions are highly 

constrained by symmetry properties and where most of the parameters can be determined experimentally or estimated 
theoretically. Therefore, for every quantity in QCD there is accessible information at very high and very low Eu- 
clidean momenta. It therefore seems like an ideal setting to use rational approximants as interpolators to fill in the 
unknown region. 

The applications of Pade approximants to the study of the strong interactions is by no means a new topic. There 
was a big surge in the 1970's, in the pre-QCD times, where multiple models of hadrons existed. With the advent of 
QCD (and the realization that the theory was asymptotically free), perturbative computations took over those hadronic 
models. Nowadays not much remains of these pre-QCD applications in reference textbooks. 

New non-perturbative techniques were developed after QCD to extract information from known low and high energy 
input. However, their motivations did not come from rational approximants and, even though essentially all of the 
methods can be recast as Pade approximants, the QCD community is only slowly acknowledging this close connection 
with Pade theory. 

In this article I will briefly review those non-perturbative standard methods, placing emphasis on their relation 
to Pade approximants. Then I will describe in some detail a set of applications to open problems where more un- 
conventional approaches based on Pade theory have been attempted and where further work in that direction seems 
promising. Let me state from the beginning that this review is not intended to be comprehensive. It is rather biased 
towards those aspects of Pade theory that I've come across in my own research. Likewise, the selection of applications 
is very personal. 



2. The large-A^ limit of the strong interactions 

An obvious possibility to gain insight into the strong interactions is to play with the parameter space of the theory. 



However, the problem is that there is no obvious free parameter other than the quark masses. G. 't Hooft [48] proposed 



in 1973 to consider the family of QCD-like theories with generic gauge group S U (N c ). In the joined limit 

N c -> oo, g 2 N c = ct. (5) 

it can be shown that a consistent non-perturbative theory of hadrons emerges. The argument is very simple. Let 
us consider, for concreteness, a generic two-point correlator with arbitrary sources. Its full result in QCD can be 
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expressed (formally at least) as a series expansion in the strong coupling constant. Whether this infinite series is 
summable is not important at this point. The important point is that the full perturbative expansion has to be dual to a 
hadronic realization in order to preserve unitarity. 

It is standard to represent the perturbative series pictorically as Feynman diagrams, some of which are shown in 
the first line of Fig.Q] The argument can be best understood by introducing the double-line notation, as shown in the 
second line of Fig.Q] where gluon lines are depicted as double lines (because they belong to the adjoint representation 
of S U{N C )) and quarks as single lines (because they belong to the fundamental representation). Double-line notation 




+ 



- 
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4^ = ^ + non — planar 

Figure 1 : Pictorical representation of a generic two-point correlator (left-hand side) in terms of its Feynman diagram perturbative series expansion 
in both conventional (first line) and double-line (second line) notations, together with its dual hadronic representation in terms of the topological 
expansion (last line). The generic sources are represented by crosses. 

thus offers an easy way to keep track of the colour flow inside a given diagram^ If confinement is assumed, then no 
open colour lines can exist. With the double-line notation it is not difficult to see that the 1 /N c expansion is actually a 
topological expansion, where planar diagrams dominate. Without solving the full set of planar diagrams, the optical 
theorem allows to connect the quark-gluon and the hadronic picture: the dashed line crossing the two-gluon correction 
in the second line of Fig.Q] gives a single colour-singlet qq state, i.e. a single meson exchange. In fact, any possible 
way of cutting any planar diagram will lead to the same conclusion [51]. Always at a qualitative level, one can show 
that the hadronic picture dual to planar diagrams consists of the exchange of an infinite number of single and stable 
resonances. In a more mathematical language, this implies that the correlator is meromorphic, consisting of an infinite 
number of single poles. Therefore for the generic two-point correlator II(g 2 ), its spectral function reads 

DO 

-lmIKt) = Yf^(t-ml), (6) 
n z — ' 

n 

where m 2 are simple poles and / 2 their associated residues. This spectral representation can be related to the full 
correlator by using Cauchy's theorem on the Argand diagram of Fig. [3] which shows the singularity structure of a 
generic two-point correlator. Therefore 

C°° dt 1 °° f 2 

n ^= —2~ t ~ Im n(0 + ^iq 1 ) = y /" 2 + Wiq 1 ) , (7) 
Jo -q 2 + tn ^ -q 2 + m 2 

where 3[{q 2 ) is a polynomial that depends on the convergence properties of II(<7 2 ) at infinity. 

Unfortunately, an analytical solution of large-A^ c QCD is still lacking. However, even being qualitative, the mero- 
morphic representation of Eq. ^ for the hadronic spectrum is currently the best available description of the hadronic 
world, and is at the base of almost all the phenomenological methods developed so far to study non-perturbative QCD. 

3. Non-perturbative methods as Pade approximants 

As already emphasized in the introduction, many of the non-perturbative methods in strong interaction physics 
developed after the advent of QCD are, sometimes in a disguised way, realizations of Pade approximants. The under- 
lying reason is that all the methods are based, to a certain degree, on meromorphic ansatze subject to some matching 
to low and/or high energy QCD inputs. In this section I will review some of the most well-established ones. 



2 group indices in SU(N C ) are normally referred to as colour indices. 
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As an illustrative example for all the methods, I will consider the following correlators: 



nj v (<7) = i j d 4 x e*'( 1 T{ V,(x) V?(0) }|0> = (<fq v - q V) iW) 

n£.(«) = ' \ d ' x eigx < 1 r ( M x ) A l(0) }| o > = («v - <?Y V ) n A ( ? ; 



(8) 



2 ) , (9) 

where V^x) = u(x)y fl d(x) and A^ix) = ii(x)j M jsd(x) are QCD currents and the factorization of the tensor structure 
follows from Lorentz and gauge symmetries. Each of the previous correlators diverges logarithmically at large values 
of q 2 , where perturbation theory can be used. However, the coefficient in front of the logarithm is the same for the 
vector and axial channel, which is a consequence of chiral symmetry being unbroken at the quark-gluon level. Thus, 
the difference Ylyv - ^Iaa = FLlr is only sensitive to non-perturbative physics and in particular to the breaking of 
chiral symmetry at the hadronic level. One normally denotes such quantities as order parameters of chiral symmetry 
breaking. Therefore, at high momenta the correlator can be parameterized entirely by the OPE as 

oo 

lim n«(^) = V^|. (10) 

»»-»(—) fri q 2 " 

It can be shown that for this particular correlator the first two terms, £2 and £4, cancel to a very good approximation. 

As a common background, all the methods to be described adopt the large-A^ f limit version of QCD described in the 
previous section. Therefore, the spectral function for YIlr can be written as 

. 00 00 

- Im U m (t) = -f 2 5(t) + J] f 2 n 6(t - m 2 Vn ) - £ f 2 n 6(t -m 2 M ), (11) 

n n 

and the correlator a^| 

„2 OO rl OO rl 

n LR (q) = — + ) — - — —-), — - 7 — 2~ ' ( 12 ) 

q 2 ^-q* + m 2 Vn ^-q 2 + m 2 An 

where the presence of the massless pion mode is a direct consequence of Goldstone's theorem. 

The differences between the different methods to be described come from the approximations made on the previous 
equation, the degree of matching to the underlying information from the fundamental QCD theory, and the different 
non-perturbative quantities one wants to compute. 

3.1. Vector meson dominance 

Vector meson dominance (VMD) l46ll rests on the hypothesis that whenever vector mesons are allowed by quantum 
numbers to contribute to a certain process, then (i) they will dominate over other particle channels and (ii) their 
lowest-lying states (the poles closer to the origin) will give the bulk of the non-perturbative effects. 

Therefore, the spectral function of Eq. (0 will be completely saturated by a single resonance on each channel, 
namely 

f 2 f2 f 2 

LR(q ) = -T + — 5 y -, J ■ (13) 

q l -q 2 + my -q 2 + m 2 A 

An additional key ingredient of VMD is that the poles and residues of the previous equation are identified with 
the physical parameters of the lowest-lying states in the hadronic spectrum, i.e. the p(770) and ai(1260) mesons. 
Therefore, in the ansatz above, my = m p = 770 Me V, m^ = m a \ = 1260 MeV, fy = f p ^ 218 MeV and = f a \ — 174 
MeV El. 



3 Due to the super-convergence of flu; at high energies, the function 3\(q~) actually vanishes. 
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3.2. Sum rules: Finite energy and Borel-Laplace sum rules 

An important problem with the VMD approach is that it requires precise knowledge of physical parameters, which 
typically is scarce, especially for the residues. The sum rule approach [47] overcomes this drawback, because it 
computes residues and/or poles in terms of a few universal OPE terms. 

The physical picture behind the sum rules is that the hadronic world, represented by the ansatze for the spectral 
function, should match QCD at high energies. Recall that in order to carry out this program one needs in general an 
infinite number of resonances in the spectral function, otherwise the logarithms of pQCD cannot be reproduced. In 
practice, and in order to be predictive, the spectral ansatz for a generic correlator Yl(q 2 ) consists of simple poles at m\ 
and a cut starting at a certain threshold Sq, i.e., 

N 

Im 11(f) = ^ f n 6(t - m 2 n ) + A6(t - s ) . (14) 

n 

Notice that despite being more predictive, this modeling of the continuum as a Heavyside-f? function departs from the 
large-A^ c version of an infinite set of Dirac-5 functions. 

There are clear advantages when one uses sum rules on order parameters. First and foremost, since perturbation 
theory cancels to all orders, one expects to be more sensitive to non-perturbative effects. Note also that in principle 
one could dispense with the continuum contribution without conflicting with the OPE. 

For ITlr, as already mentioned, not only perturbation theory cancels, but also the terms q with d < 6 in the OPE. 
This means that the matching equations read 

N A N v 

ft = - £° *imn„(o = 2>a 2 «-Z^» + ^ 2 s o 

n n 
N A N v 

ft = ~ $™ dttlmn LR (t) = Yjti" m2 *"-Yj f V>< m Vn = 



N A N v 

2 2 7 -2 
Vn m Vn 



where £2/ are the OPE coefficients defined in Eq. JlOt . In the sum rule approach the poles are taken as input from 
the experimental masses of the lowest-lying states in the spectrum. With ITlr, experience indicates that the matching 
equations above can be satisfied with the minimal hadronic content. However, in certain cases this might turn out to be 
a bad approximation and contributions from higher order poles can be important. Instead of adding those contributions 
explicitly, a common strategy is to use Borel sum rules. 
With Borel sum rules, one starts from 



X 



dte-' T - Imn LS (?) = -— (h e- c > T n LR (q 2 ) , (15) 

71 lm J\q 2 \=s 

which again follows from Cauchy's theorem applied to Fig. [3] where so is the radius of the circle. On the left-hand 
side one plugs in the hadronic ansatz while on the right-hand side one uses the OPE0 In general one will find 

n n j 

Notice that higher values of the Borel parameter t suppress exponentially the contributions from heavier states on the 
hadronic side, but at the same time enhance the weight of higher dimension operators in the OPE side. In practice one 
should look for some stability (i.e. r-independent window), in the Borel plane. 



4 Strictly speaking, plugging the OPE is inconsistent with Cauchy's theorem. In general, on the right-hand side one should also include terms 
that account for quark-hadron duality violation. I will not delve into this issue and instead refer the reader to [9], where a systematic study of such 
effects on the fltR correlator is performed. 



P2 Xn LR (Q 




Figure 2: Compai'ison of the MHA ansatz (solid red line) with experiment (dotted green line) on IIlr together with the predictions of other methods 
without high energy matching. Figure borrowed from Ref. I 44ll . 



3.3. The minimal hadronic approximation 

The minimal hadronic approximation (MHA) 11391 12711 was initially devised to compute non-perturbative quantities 
that can be expressed as integrals over the Euclidean regime of correlators. One such instance is the pion electromag- 
netic mass difference, given by [ 14] 



tsm„ 



3a r 

Jo 



dQ 2 Q 2 Tl LR (Q 2 ) 



(Q 2 



-q 2 )- 



(17) 



The MHA assumes the correlator to be meromorphic, with the minimal number of poles to reproduce the OPE behav- 
ior known from QCD. It rests on three key points: 

• the importance of matching to both low and high energies, i.e., to the Laurent (OPE) and Taylor (chiral) ex- 
pansions of the correlator under study. There are different reasons to prefer this 2-point matching to a simple 
1 -point one: first, the more information one can implement in the ansatz the better, and the first terms in both the 
Taylor and Laurent expansions are known with reasonable accuracy. Second, requiring matching on both high 
and low energies improves the stability of the result, especially since one is interested not in the local behavior 
of YIir but on the area below it. See for instance Fig.|2]borrowed from Ref. [44]. 



• The poles of the meromorphic function are identified with physical particle masses, while the residues are 
determined from the matching equations. 

• The method is restricted to correlators that are order parameters, i.e., those for which the contribution from 
pQCD identically vanishes. In other words, the meromorphic ansatz is only compatible with ultraviolet finite 
correlators, namely those that behave like inverse powers of q 2 ^ Correlators with Euclidean logarithms are 
therefore beyond the scope of the method. In practice however this is not a severe restriction on the number of 
observables to compute. 

For further details and different applications of the MHA to electroweak physics I refer the interested reader to 

Ref s. EH [H m [M Si . 

From the previous discussion it is evident that both sum rules and the MHA are straightforward applications of 
Pade-type approximants. FESR are 1 -point Pade-type approximants, BSR are 1 -point Borel-Pade-type approximants 



5 Even in this case the method has its limitations. For instance, it has been realized [4] that a meromorphic ansatz with a finite number of poles 
cannot fulfill all the high energy constraints, especially those coming from quark-counting rules. 
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while the MHA is an application of 2-point Pade-type approximants. Nonetheless, let me emphasize once again 
that all the previous methods did originate without prior knowledge of Pade theory. Actually, in all of them, the 
parameters appearing in the ansatz (poles and residues) are considered as physical parameters, something unjustified 
from the Pade theory perspective, where the ansatz is seen as a mere interpolator. I am fully convinced that these 
prejudices behind the particle interpretation of the ansatz are precisely what has prevented to move from Pade-type to 
regular Pade approximants in non-perturbative QCD. After all, it is well known that Pade poles are not restricted to 
be real-positive. In any case, the fact that the previous approaches turned out to be realizations of Pade approximants 
is not only remarkable per se, but at the same time makes Pade approximants compelling as natural tools to address 
non-perturbative QCD problems. 



4. Some applications 

In this section I will describe in some detail three different open problems in non-perturbative QCD that cannot 
be addressed with the techniques reviewed in the last section. However, the potential of Pade approximants to make 
progress in each of them is far from being exhausted. 

4.1. Pade approximants and the hadronic spectrum 

At the end of the 1970's there was an attempt to determine the spectrum of QCD (or at least the one of large-A^ c 
QCD) using Pade approximants. To the best of my knowledge, this is one of the first applications of Pade theory to 
non-perturbative particle physics after the advent of the QCD Lagrangian. In the following I will briefly outline the 
original derivation of Ref . B36I1 . 

The Pade approximant to Ylyiq 2 ) (cf. Eq. l[8]l) around a point is given by the following equation: 



UNiq ) 



where by construction U}y' M \q 2 ) reproduces the first N+M+ 1 derivatives of Yly(q 2 ) around the point. In the following 
I will consider the diagonal sequence [N, N]. Hence it follows that 



d" 



Tl v (q 2 )Q N (q 2 )-P N (q 2 ) 



= 0, n = 0,...,2N. (19) 

cf=-fi 2 



d(q 2 )" 

By construction, cancels after the first N equations, and applying Cauchy's theorem one ends up with 

dt 



f 

Jo 



. , ■ lmTl v (t)Q N (t) = 0, n = N + 1, 2N , (20) 

Jo (t + p z ) n+[ 

from which one can determine the denominator Q^. So far no approximation has been made. In order to proceed 
further and solve the previous equation some information has to be provided on ImTlv(t). Back in 1977 very little was 
known about non-perturbative expansions (the operator product expansion is from 1978), so the natural starting point 
was the perturbative expansion, whose leading term reads 

2,_ N c ^J-q 2 



iim n v (^) = -prolog Mr ' (21) 

9 2 -»(-°°) 127T 2 \ p z j 

This selects q 2 - -p? as a natural expansion point, depicted in Fig. [3] With the previous input, Imny(f) is simply a 
constant and Eq. ( f20b turns out to have a simple solution, namely 

<W) = *t (-N, -N;i;- q ^) = + pY pT . (22) 

where 2F\{a,b;c\d) is Gauss' hypergeometric function and P^ )0) stands for the (0, 0)-Jacobi polynomials, i.e., the 
Legendre polynomials. Plugging the solution back in Eq. ( fT9l one can determine 'P^iq 2 '). The result for the Pade 
approximant eventually reads [49] 




Figure 3: Singularity structure of two-point correlators TL(q 2 ) in the complex q 2 plane with the contour chosen to apply Cauchy's theorem. The 
circle is defined by \q 2 \ = so, which will be eventually sent to infinity. The singularities on the right half-plane are the poles ofYl(q 2 ), sitting at m 2 . 
The point on the left half-plane is the pole with multiplicity (n + 1) explicitly shown in the denominator of Eq. (20). 



The problem therefore has reduced to computing the Pade approximant to the natural logarithm, a result originally 
found by Jacobi, Rouche and Gauss back in the 19th century. 

Convergence of the Pade approximant means that the original function is recovered when one lets N — > oo. In the 
case of the logarithm it is well known that the set of Pade poles becomes dense over the physical axis eventually 
mimicking the logarithmic cut. However, the main motivation behind Ref. ll36ll was to extract a set of poles from a 
continuum and a Pade clearly does not serve this purpose. The following correlated limit was taken instead, 

q 2 « //, N oo, — = ct. , (24) 

leading to the following difference in the resulting functions: 

[N -> oo} , (25) 

\N -> oo, q 2 « //, — = ct.) . (26) 

Notice that even after the continuum limit is taken, in the second line there are a set of poles sitting at the zeroes of 
the Bessel Jo function (see Ref. [7] for a discussion on the physical interpretation of this last term). Those poles were 
claimed to give the spectrum of Tlyiq 2 ), in other words, the spectrum of vector mesons. 

There are a set of criticisms one can make to this approach, the most obvious being that the whole procedure is 
not a Pade approximant, contrary to what it was initially advertised. On the physics side, some comments are in 
order too. First and foremost, I find it difficult to justify that starting from the perturbative quark-gluon logarithm 
one can infer the spectrum of hadrons, a purely non-perturbative entity. Second, but closely related, constructing 
the Pade around an asymptotically large point fi 2 is problematic: the expansion at large momenta is an asymptotic 
expansion of the original function, hence it does not uniquely define the original function. In other words, the resulting 
meromorphization of the function is not unique. In fact, one can easily show that different spectra can reproduce at 
high energies the same partonic logarithm |7j]. This would be avoided if the Pade is constructed around a finite point, 
for instance q 2 — l38ll . where data exists and the chiral expansion is well defined. 

Although the solution proposed by Migdal turns out to be wrong for several reasons, the problem is very interesting 
and its solution not yet settled. I honestly believe that a proper application of Pade approximants can improve the 
current situation. 

For instance, one could start by feeding Eq. d20b with OPE condensates, which would certainly incorporate non- 
perturbative physics into the problem. In the same spirit, one could argue that some (Borel) resummation of the 
perturbative series might also help: to the extent that the asymptotic character of the perturbative expansion is thought 
to be a consequence of the existence of non-trivial infrared physics, there is definitely non-perturbative information 



n N v (q 2 ) 



iW) 



4 N e , 

"w log 



Q N (q 2 ) 



4 N c 
~3 (47r) 2 



q 2 Y (qy m ) 
log — - n— 

LL l Jo (qy m ) 
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lurking behind pQCD. This is at the base of the theory of renormalons (see the contribution of Prof. J. Fischer in this 
conference). 

While I cannot exclude that this approach can yield some benefits, I foresee two potential difficulties: first, the 
perturbative expansion is plagued with logarithms, which cannot be used as input in a conventional Pade approximant; 
and second, even if the previous difficulty is circumvented, one has to make sure that the Borel resummation captures 
a substantial amount of non-perturbative effects. 



What seems more feasible to me is to start from an expansion around the origin, in the spirit of 13811 . Then no 



logarithms are present, and the non-perturbative physics can be easily implemented from the following moments: 

-lmn y (0=^, (27) 
o t 

where £ n are the coefficients of the MacLaurin expansion of Ily, which can be extracted from experiment. Addition- 
ally, it can be shown that Ily is a function of Stieltjes type [38']- Therefore, the Pade poles all lie on the physical axis 
and convergence of the Pade approximant (at least for finite q 2 and out of the physical axis) is guaranteed. 

I do not know how far this approach can go. Certainly there are theorems like Koenig's theorem 12911 to approximate 
the poles of a function. The problem in QCD is that the number of poles is infinite, and I do not know of any result 
that applies in this case. 

4.2. Pade approximants and the gauge/string duality 

Quite recently there have been attempts to infer the spectrum of QCD from theories defined in (4 + d) dimensions. 
The idea that our space-time, i.e., a 4-dimensional manifold with metric rj^ = diag(l, — 1, -1, -1), and the quantum 
fields defined in it could be dynamically generated from a higher dimensional manifold is quite old. More recently, 
this kind of constructions have been suggested by superstring theories, where the cancellation of anomalies requires 
theories to be constructed in ten or eleven dimensions. Since in our everyday life we only perceive 4 dimensions (3 
space-like and 1 time-like), the remaining extra dimensions have to be dynamically compactified. What has made 
those ideas more compelling are the recent developments following the so-called AdS/CFT correspondence (see for 
instance Ref. yfl and references therein), that gives a precise and well-defined prescription to connect a 5-dimensional 
manifold endowed with Anti-de Sitter (AdS) metric 

ds 2 = g MN dx M dx N = Xr{-dy 2 + Tj^dxf 'dx v ) , (28) 

y 



to a 4-dimensional submanifold, through the following identity [53]: 



Ssd[4>o\ = J d 5 x£ 5D (<p Q ) = J d 4 x(£ 4D + <p a Oi) = S 4 d + sources , (29) 

where 0o are the solutions of the 5-dimensional equations of motion projected on the 4-dimensional manifold and (9,- 
are operators in the 4-dimensional theory. Therefore, the correspondence claims that the 5-dimensional field solutions 
on the 4-dimensional submanifold are the sources of the correlators in the 4-dimensional theory. Even in these rather 
exotic scenarios Pade approximants might help. 
Let us concentrate, as before, on vector mesons. We start from the massless Yang-Mills action in five dimensions: 

S = ^| J d 4 x £' dy^Tr[F MN F MN ] , (30) 

where the 4-dimensional boundary branes are located at y — e — > (high-energy or ultraviolet brane) and y = y m 
(low-energy or infrared brane). The setting is depicted in Fig. [4] Without loss of generality one can set V5 = 0. Then 
the equation of motion for the field reads 



1 



^--^-□1^=0, (31) 
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Figure 4: Geometric setting used in Section l4~2l The left-hand side manifold is the UV brane, where the gauge theory is defined. The right-hand 
side brane (IR brane) defines a particular model of confinement. The space in between, also called the bulk, is a 5-dimensional manifold with 
Anti-de Sitter metric. 



which, subject to the following boundary conditions 

dyVffym) = 0; 



results in [42] 



(32) 



(33) 



e Y (qy m )Ji(qe) - J (qy m )Yi(qe) 
Following the AdS/CFT prescription, the previous solution can be plugged back into the action, leaving only a bound 
ary term 



2g 2 J 



d 4 x -VudyV^ 
y ' 1 



(34) 



Notice that vf \ which is left unspecified, is the value of the field V u at the UV brane, and therefore also the source of 



correlators in the 4-dimensional theory. Therefore, Yly is given by [ 17 

id id „ 



5^ SVf* 



(q"q v - tfsT) 



l l 



dyV 



(35) 



The solution of the previous equation is exactly the one given by Eq. (l26b in the previous section, with the identification 

2N_ 



y« 



(36) 



In other words, Migdal's prescription and the holographic model described before turn out to describe the same 4- 
dimensional theory. This has led some works ll 811 to erroneously conclude that the Pade approximant to the spectrum 
of QCD arises naturally in 5-dimensional models with AdS metric. However, after the discussion of the last section, 
this statement cannot be right, as first pointed out in 01. 

The reasons why both approaches yield the same results have been discussed in detail in Ill9h , Roughly speaking, it 
has to be attributed to three common ingredients, which are implemented in very different ways: 
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• conformal invariance. The Lagrangian of QCD in the massless limit is scale invariant. Thus, conformal 
invariance is implemented by hand in Migdal's construction by choosing the leading perturbative logarithm as 
input. In the holographic model the presence of the logarithm is just an automatic consequence of conformal 
invariance of the AdS metric. 

• infinite spectrum. Through meromorphization in Migdal's construction, and by compactification in the holo- 
graphic model. Actually, the 5-dimensional equations of motion for the Vf, field can be cast as a quantum 
mechanical eigenvalue problem. 

• confinement scale. Introduced in Migdal's construction through a sophisticated continuum limit, and set by the 
infrared brane in the holographic model. 

Albeit at present both approaches are rather ad hoc mechanisms to mimic confinement, the similarities between them 
are extremely interesting. Applications of the AdS/CFT correspondence to QCD require a departure of the conformal 
limit in order to take confinement into account. If one knew how to implement this in a solid way, one would have 
valuable insight on how the string dual of QCD looks like, one of the holy grails in the string community. If progress 
can be made on the 4-dimensional side with Pade approximants (perhaps along the lines discussed at the end of the 
last section), then the gauge/string duality could be used to make contact with string theory constructions. 

4.3. Pade approximants in Euclidean space: relations between low and high energy parameters 

I will finally report on a different application of Pade approximants, namely the prediction of high energies from 
known low energies. For the sake of illustration, I will concentrate on the Wir correlator, with generic meromorphic 
ansatze of the form 

rl M' fl Na fl 

V. m (q ) = — + > — 5 — ~ ), — 5 — > (37) 

a z ^— ' -q 1 + mi -q 1 + mj 

2 n . " vn n " An 

The discussion will follow closely Ref. fl. 

I already mentioned that at high energies, the first two terms of the OPE vanish and therefore the correlator converges 
like q~ 6 , meaning that £2 = and £4 = 0. A determination of the leading fa and subleading coefficients of the 
OPE from experiment turns out to be problematic, and different analyses find serious discrepancies (see the first two 
columns of Table 1, where existing phenomenological determinations are summarized). The only solid aid from 



theoretical considerations comes from a theorem by Witten [52], which can be cast as the inequality 



q 2 n LR (q 2 ) > 0, -00 < q 2 < . (38) 

The previous equation in particular implies that ^ > 0, but nothing is known about Actually, notice from Table 1 
that even its sign is disputed. 

However, it turns out that rather accurate experimental data exists on Imn^f) at low energies and in principle 
the chiral parameters (the MacLaurin coefficients) of TIlr can be determined with reliability. Taking q 2 TliR(q 2 ) as 
given by Eq. (l37l i as a Pade approximant, relations between high and low energy parameters can be established. The 
minimal version of Eq. (|37T i consists of just one vector and one axial contribution, subject to the following constraints: 



x 

Jo 



dtlmIl LR (t) - 


- /a fv + fn ~ 


0, 


dttlmTl LR (t) - 


= fWk - fWv 


= 0, 


dt 

-Imn„(?) = 


fl fl 

Jv Ja 
ml m\ ' 




dt 

T Imn LS (f) = 

t L 


f i f i 

Jy J A 
4 4 ' 

my m\ 





? 4 = I ^Imn is (f) = ^t-^t, (39) 
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where the first two equations are the high-energy constraints and the last two are the low energy matching equations 
coming from the MacLaurin expansion (chiral expansion) of the correlator, defined as 

f 2 

limIl LR (q 2 ) = ±§ + Y C 2j q 2j - (40) 
q 2 ^o q z t-r 1 



j 



The parameters fy, /a, my and rtiA are therefore determined from Eqs. d39l l as functions of f„, t,\ and £,%. The expres- 
sions for ^6 and £g turn out to be rather simple JH] : 

£i - (aft; 



There are a set of consequences that can be readily inferred from the previous equations: 

• Since it is experimentally established that > 0, Eq. (l42l immediately implies that fg > 0. Therefore, the 
Pade approximant we are considering here favours the set of phenomenological analyses sitting on the first half 
of Table 1. 



Combining Eqs. (PTTT i and ( l42l one can get the relation 



6. = , (43) 

Jn 



which turns out to be fulfilled by the same first half of Table 1 to a remarkable degree of accuracy, as shown in 
Fig. [5] Therefore, even though the different phenomenological determinations differ on the values of ^ and £g, 
it seems as if Eq. (l43l is a universal constraint on those determinations with £g > 0. 



6<It-I ' ;44s 



Witten's inequality on Eq. ( PiTT i sets an upper bound on 

\.fn 

which can be used as a consistency check of the full approach. Plugging typical numbers for £i and f„, one gets 
^ 2 < 0.1 GeV" 2 , which agrees well with the typical & 0.08 GeV ' 2 . 

The exercise considered before is just the first step of an iterative process, where one should consider the sequence of 
Pade approximants P% +2 subject to an increasing number of low energy input in the form of the MacLaurin coefficients 
%2j- This will give rise to an associated sequence of predictions for ^ and £s- What is remarkable is the non-trivial 
agreement with the first half of Table 1, which suggests that already P Q 2 might yield values for and £s close to the 
real ones. Should this be true, then we would have a criteria that strongly favours > 0. 

This kind of strategy is especially suited for problems where only Euclidean quantities, i.e., those defined or involv- 
ing only the Euclidean half-plane, are requested. Notice however that this approach differs from the MHA introduced 
in Section 13.31 The MHA is an exercise in Pade interpolation, while the approach discussed here is an exercise in 
Pade extrapolation. In fact, extrapolation from q 2 — to arbitrary large values of q 2 ^ Notice also that I use Pade 
approximants instead of Pade-type ones. 



Notice that the problem falls out of the scope of Pommerenke's theorem 1 43] because the OPE is defined at q —> oo, i.e., not in a compact 
subset around the origin. 
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4_Q 71 4- f) OA 

+ 7. / 1 It U.7U 
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+fi + 6 


4_7 6 + 1 5 


+5.2 +1.0 


Cirigliano et al. on ALEPH [13 ] 


+4.45 ± 0.70 


-6.16 ± 3.11 


+2 86 + 90 


Cirigliano et al. on OPAL [13 ] 


+5.43 + 0.76 


-1.35 + 3.47 


+4.3 ± 1.2 


Bijnens et al. on ALEPH [5J 


+3 4 +2A 

" -Z0 


_i 4 4+10.4 


+ 1.7 + 2.4 


Bijnens ef a/, on OPAL [5] 


+4.0 ± 2.0 


-10.4-1;° 


+2.3 + 2.3 


Latorre et al. f451 


+4.0 ± 2.0 




+2.3 ± 2.3 


Almasy et al. [3j 


+3 2 +L6 


-17 +25 


+ 1.5 + 1.5 


Ref. [8] 


+7.6 ± 0.4 


+8.3 ± 1.0 





Table 1: Values for the dimension-six and dimension-eight OPE condensates (in units 10~ 3 GeV 6 and 10~ 3 GeV 8 , respectively) reported using 
different phenomenological techniques. In the last column I list the would-be value for the dimension-eight condensate if Eq. 1431 were used, taking 
as inputs typical values for and f n and the different values of 



5. Conclusions 

I have presented a brief overview of commonly used methods to compute non-perturbative quantities in QCD and 
their close relationship with Pade approximants to meromorphic functions. All those methods were originally devel- 
oped without prior knowledge of the developments in Pade theory. Only recently this connection has been acknowl- 
edged and its far-reaching consequences are only starting to be fully developed J3~ll l8ll. It is quite remarkable how 
concepts like meromorphization come out rather naturally in QCD if one embraces the large-A^. limit; quark-hadron 
duality guarantees the matching between pQCD (and the OPE) and the hadronic spectrum; or order parameters of chi- 
ral symmetry breaking avoid the presence of high energy logarithms. The fact that the different physical approaches 
ended up converging to well-established methods in numerical analysis reveals Pade approximants as extremely useful 
tools to explore the non-perturbative aspects of QCD. 

Certain problems, like the ones I presented in Sections [4. 1114. 2l and l4. 31 fall out of the scope of such techniques and 
require new methods. I am convinced that strategies based on Pade approximants can lead us far, but attempts have 
been scarce so far. 

In writin g th is article I had to omit quite a lot of material, such as the applications of Pade theory to heavy quark 



physics (6i[12j,|24j,|33|], the interesting works on light quarks of Refs. [32, 34], Pade unitarization Il35ll or QCD with 
finite temperature and chemical potential 13011 . Hopefully the contributions of S. Peris and J. J. Sanz-Cillero in this 
conference can help fill some of this gaps. 
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